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Abstract 

Relying on a few lowest order perturbative calculations of anomalous dimensions of gauge 
invariant operators built from holomorphic scalar fields and an arbitrary number of covariant 
derivatives in maximally supersymmetric gauge theory, we propose an all-loop generalization of 
the Baxter equation which determines their spectrum. The equation does not take into account 
wrapping effects and is thus asymptotic in character. We develop an asymptotic expansion of the 
deformed Baxter equation for large values of the conformal spin and derive an integral equation 
for the cusp anomalous dimension. 



1 Introduction 



Four- dimensional nonabelian gauge theories were found to possess integrable structures. Initially 
they had been discovered in evolution equations for scattering amplitudes in multi-color limit of 
high-energy QCD [H [2]. Subsequently it was demonstrated that hidden symmetries also arise 
in the field-theoretical dilatation operator whose eigenvalues determine anomalous dimensions 
of gauge invariant composite operators of elementary fields in underlying models. Integrability 
was revealed in one-loop anomalous dimensions of twist-L maximal-helicity quasipartonic Wilson 
operators in QCD |3j by identifying the former with eigenenergies of the L— site XXX Heisenberg 
spin chain |3]. The magnet turns out to be noncompact, for the spin operators acting on its 
sites transform in the infinite-dimensional representation of the coUinear subgroup SL{2, M) of 
the conformal group SO {4, 2). Since the one-loop phenomenon is spawned by gluons, invariably 
present in Yang-Mills theories, — supersymmetric or not, — they all necessarily exhibit the same, 
universal integrable structures. The differences arise due to distinct particle contents of the 
models: while only holomorphic sectors are integrable in QCD and its nearest supersymmetric 
A/" = 1, 2 siblings [5], the maximal supersymmetry of the A/" = 4 super- Yang-Mills theory extends 
integrability to all operators PIZIIHIE]. Recent perturbative studies build up a growing amount 
of direct evidence that integrability persists in certain closed compact [3, [8] and noncompact 
P, [ini [m [121 Us] subsectors of gauge theories even in higher orders of perturbation theory. 
Thus, while ruled out for gauge theories with Af < 4 supercharges, it is plausible that the max- 
imally supersymmetric Yang-Mills theory is completely integrable. An additional confirmation 
for this conjecture comes from studies of multi-loop multi-leg scattering amplitudes which display 
intriguing iterative structures [HI [I5]. These arguments suggest that the spectrum of all-loop 
anomalous dimensions in AA = 4 SYM theory is determined by a putative long-range integrable 
spin chain with the dilatation operator being its Hamiltonian. 

In this note we probe the underlying integrable long-range magnet by proposing its multi-loop 
perturbative structure within the framework of the Baxter Q— operator [16]. This approach is 
based on the existence of an operator Q{u) depending on a spectral parameter u and acting on 
the Hilbert space of the magnet. For different values of u it forms a family of mutually commuting 
operators, simultaneously commuting with the spin-chain Hamiltonian as well. Although in the 
present circumstances, the formalism is equivalent to the Bethe Ansatz approach, it possesses 
certain advantages. First, the eigenvalue Q{u) of the Baxter operator Q{u) determines the single- 
particle wave function of the chain in the representation of separated variables [T7] . Second, the 
equation for the Q— operator, — known as the Baxter equation, — is polynomial, to be contrasted 
with a set of coupled transcendental Bethe equations. Third, it allows for a straightforward 
asymptotic analysis when quantum numbers of the chain are large as will be demonstrated 
below. 

Currently we restrict our consideration to the closed noncompact SL{2) sector of the gauge 
[TSl [TU] theory which is spanned by single-trace maximal i?— charge Wilson operators built from 
the holomorphic scalar fields X = 01 + i(j)2 and covariant derivatives, 

0n,n....n,(O) =tr{(^D+rX(0)(zD+)"^X(0)...(^Z}+rX(0)} . (1.1) 

Here Dj^ = D^n^ is projected on the light cone with a null vector ra^, = 0, in order to factor 
out the maximal Lorentz-spin component from the operator in question. These Wilson operators 
mix with each other under renormalization group evolution and acquire anomalous dimensions 



at all orders of perturbative series in coupling constantQ 

oo 

7(^) = E^''^^^"^- (1-2) 

n=l 

We find it convenient to use the expansion parameter g related to the 't Hooft coupling constant 
A via 

g = ^2\ = ^^. (1.3) 

The anomalous dimension "^{g) depends on parameters characterizing the operator: its twist L, 
determined by the number of X— fields, and its Lorentz spin N = rii + n2 + ■ ■ ■ + n^. Within 
the method of the Baxter Q— operator, the eigenspectrum of one-loop anomalous dimensions 7'°' 
and the corresponding quasimomentum are determined by the leading order Baxter function 

7'°> = f [lnQ'°)(f)]' - i [lng'«)(-f)]' , = lnQ(°)(|) - lng""(-f) . (1.4) 

Since the Baxter function Q^^\u) is related to the eigenfunction of the mixing matrix, it cor- 
responds to a multiplicatively renormalizable Wilson operator and thus has to be polynomial 
in u of order (^'"'(m) = {u — uf^){u — Mj"') . . .{u — u^^). The zeros of this polynomial are 
determined by the Bethe roots u^^'' which take only real values for the noncompact SL[2, M) spin 
chain [2l|20]. The function Q^°\u) obeys the finite-difference Baxter equation [16] 

u';_Q''\u + u'lQ^'^u - z) = t^'\u)Q^'\u) , (1.5) 

where the spectral parameter in the dressing factors is shifted by the conformal spin s = | 
of the scalar field X , u± = u ± ^ and t'°'(u) is an order— L polynomial in u depending on the 
integrals of motion. 



2 Three-loop Baxter equation 

Explicit perturbative calculations [9l [101 [12] of two-loop corrections to the anomalous dimensions 
of the scalar operators (11.11) exhibit double degeneracy of energy levels with zero quasimomentum. 
This hints at the existence of nontrivial odd-parity conserved charges and thus persistence of 
integrability at higher orders of perturbation theory. 

Beyond one loop, the formalism of the Baxter operator gets modified accordingly. The Bethe 
roots acquire corrections in coupling constant to all orders of perturbation theory, 

00 

Un{9) = J29"V, (2.1) 

and obey deformed Bethe Ansatz equations [21]. The reality of Bethe roots Uk{g) have to be 
preserved to all orders since the eigenvalue Q{u) of Q{u) is a wave function of the chain with 

^ Their complete two- loop planar mixing matrix has been recently computed in Ref. [13j . 
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the number of its nodes on the real m— axis coinciding with the spin N of the operator. The 
polynomial 

N 

Q{u) = l[{u-un{g)), (2.2) 

n=l 

fulfills these properties and is real Q*{u) = Q{u*) for u* = u. In Ref. [13] we found from available 
two- [22], [23l m [To] and three- loop [MJ |25l [26] diagrammatic calculations of anomalous dimensions 
that the Baxter equation possesses the form 

x^e"+(^+)g(M + i)+ x^e''-^''-^Q{u - i) = t{u)Q{u) , (2.3) 

with the dressing factors depending on the renormalized spectral parameter [27] 

x[u] = ^ + - , x± = x[u±] . (2.4) 
The multi-loop transfer matrij|§ 

tiu) = 2m^ + qi{g)u'^-' + q2{g)u''-^ + ■■■ + q^g) (2.5) 

acquires the "missing" term ~ u^~^ at (7^— order, i.e., ^'1(5') ~ 0{g'^), while the rest of the charges 
start from 0{g^), qk{g) = ql'^ + C'((?^). The additional dressing factors a± obey the complex 
conjugation condition ((j_(x_))* = a^{x^) for SJm-u = and encode the renormalization of the 
noncompact charges qk{g) at higher orders. An analysis yielded the following result to three- loop 
order [13] 

^±{x) = [lnQ(±f )]' - { [lnQ(±f )]" + X [lng(±f )]'"} + 0{g') . (2.6) 

While the anomalous dimension is expressed order-by-order in coupling constant g in terms of 
the solution to Eq. (12.31) as [I3] 

1(9) = ^{y [l^QW]' + ^ [lnQ(«)]'" + ^ [\nQ{u)f^ + 0{g')y . (2.7) 

The anomalous dimensions found using these equations reproduced exactly available perturbative 
predictions. One can demonstrate that the condition of the pole-free transfer matrix at Bethe 
roots Un{g), t{un) = immediately produces the three-loop Bethe Ansatz of Ref. [28] . 

3 Multi-loop conjecture 

The above representation (12.61) of the dressing factors a± can be brought to a very suggestive 
form. Namely, a quick inspection allows one to rewrite these terms as an expansion in terms of 
the Chebyshev polynomials of the second kind f/^, 

= f / * [lnO(±i - gt)]' (-A) Ml . (3.1) 

"^^1 n=0 

^Note that with this transfer matrix the resuhing Baxter equation breaks down aheady at order 0{g^"') in 
couphng constant with n — L. It turns out that one can correctly incorporate order n = L corrections by replacing 
the leading term in t{u) with the following combination 2u^ ^ + — (|)^ — (— f)^- 
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with nmax = 2, valid to 0{g^) in the approximation of Eq. (\2.Q\\ . Having this representation at 
our disposal, we may naturally extend the first few terms of the available perturbative series to 
all orders in coupling g, by sending nmax oo. Using the summation theorem for Chebyshev 
polynomials, one can sum the infinite series up into the function — (arccot^^0)/vT^^T^ ^nd, 
upon a variable transformation, write a± in the form (with z = 1 — z) 



[\ dt 



a±{x) = —- — / dz 



2vrx Jo J_i v^r^ 



Ing ( ±^-g^t + z^ 



^ (3,2) 

7-1 v^r^ u + gt ^ ^ 

Here we integrated by parts in the second line in order to separate the components 6± of the 
spin-chain quasimomentum 9 = 9^ — 9_, 

r-gl..Q(±|-.i) 

Notice that 9 reduces to the one-loop expression ( ll.4p upon setting (7 = 0. While the condition 
t{un) = yields the all-order Bethe Ansatz equations suggested in Ref. [2Tj . 

The conjectured multi-loop Baxter equation fl2.3p with fl3.2p and the known pattern of renor- 
malization of the conformal spin in field theories can be used to determine the all-loop analytic 
expression for the anomalous dimensions in terms of the Baxter function. To this end, recall 
that the conformal spin of Wilson operators J'°' = N + defining the quadratic Casimir 

= —J^°^{J^°^ — 1) — gets additive renormalization by the anomalous dimensions '-/{g) of 
composite Wilson operators at higher orders in coupling, J'"' J = N+^L+^j^g). This conclu- 
sion arises from considerations of conformal Ward identities for Green functions with conformal 
operator insertion [291 ■ Then a short inspection of the Baxter equation (12.31) with the dressing 
factors a± in the form (13. ip demonstrates that the first term in the series of (t±{x) = i'^±{g) / x+. . . 
induces the shift of the conformal spin, 

j(o) = N + \L^J = N+\L + \ (7+((7) - i-{g)) . (3.4) 

Consequently, we may naturally identify the addendum with the anomalous dimensions of a 
multiplicatively renormalizable composite operators, "^{g) = l+ig) — l-{g)- Making use of the 
explicit form of the dressing factors a±, we find the integral representation of '~f{g) in terms of 
the solution to the Baxter equation, 

lig) = I- f dt VT^ [\nQ {^,-gt)-\nQ (-f - gt)]' . (3.5) 

The Taylor expansion shows that the lowest three orders in g^ coincide with Eq. (12. 7p . 

The Baxter equation (12.30 can be solved analytically order-by-order in coupling constant for 
specific values of L and A^, e.g., for L = 4, = 2 eigenvalue with zero quasimomentum read^. 



1(9) = '-^g' - '^g' + '-^^^g' - ^^^^f^/ + 0{g'') . (3.6) 



•^This anomalous dimension, when related to Berenstein-Maldacena-Nastase operators [30], agrees with previ- 
ous one-, two- and five-loop analyses of Refs. [31], [7] and [27], respectively. 
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However, it has a limited range of applicability being asymptotic in character: it allows to find the 
anomalous dimensions up to order 0{g'^'"') only for operators of length L > n. This restriction 
arises from the breaking of its polynomiality above a boundary value of n, i.e., for L < n. 
Analogous limitations apply to the Bethe Ansatz equations of Ref. [27]. A generic dependence of 
7 ((7) on the parameters L and N is not known however and below we will develop an asymptotic 
scheme to find it in the large spin limit. 



4 Asymptotic expansion 

The large— behavior of anomalous dimension is of special interest in its own right since it 
governs the Sudakov asymptotics of scattering amphtudes [321 133], and in light of gauge/string 
duality, for it can be compared (at strong coupling) to energies of quasiclassical strings [Ml [351 
l36l [371 [38] . Recall at first that the anomalous dimensions of twist— L operators occupy a band 
of width L — 2, with the upper and lower boundaries scaling like [H [T8] 

7iower(fi') = 2rcusp(fi') In A^ , 7upper(5') = LTcnsp{g) In A^ , (4.1) 

and the coefficient rcusp(5') being the cusp anomalous dimension [391 [IQ], known to one- [39], two- 
[iQl [221 [231 [18] ^-nd three- loop orders [251 [261 [15] . The minimal anomalous dimension 7iower(fi') of 
high-twist operators develops the asymptotic behavior identical to the one of twist-two operators 
[371 [T2]. Since the single- logarithmic regime is realized for Le^ -C A^ with L, N —!■ 00 [37], this 
allows one to evade the limitation of the asymptotic character of the Baxter equation and to 
derive an all- loop equation for the cusp anomaly Fcusp- 

Notice that although we have to solve the problem with large quantum numbers, we cannot 
apply traditional WKB expansion for Q{u) (see, e.g., Ref. [4]) since the latter is valid for the 
spectral parameter which scales as m ~ A^^ while the energy is determined by the Baxter function 
Q{u) evaluated at the argument u = ±^ — gt which behaves as m ~ A^°. Therefore, we have to 
resort to other techniques. To this end, we will use in the following the approach developed in 
Refs. [^ [37] for one-loop anomalous dimensions and which, as we will see momentarily, is easily 
generalizable beyond leading order of perturbation theory. 



4.1 One- loop Baxter equation 

Let us briefiy review the formalism of Refs. [HI [37] applied to the one-loop Baxter equation (11. 5p . 
Though we are interested only in the lowest energy curve, at the beginning we will be general 
enough to discuss subleading trajectories as well in order to point out approximations which have 
to be imposed to separate the lowest anomalous dimension only. In the regime in question, the 
conserved charges are large ~ A^'^ and, therefore, the transfer matrix is large |t'-°'(M)| ^ 1. 
Introducing a new function 



we can rewrite the Baxter equation in the form 

0(0) (iT- i) 



n^0<°)(n) + -— — — = t'°'(n). (4.3) 
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The solution to it is based upon different scaling behavior of the right- and left-hand sides with 
N. For the spectral parameter u ~ A^*^, the solution is given by an infinite fraction. Keeping the 
leading terms only we come to two difference equations 

«+<5+ {u + t)= {u)Ql^ (u) , u'lQ^'} {u-i)= {u)Q'°^ (u) . (4.4) 

The additive corrections to their right-hand sides go as where is a conserved charge 

which scales with the maximal power of A^. For cyclically symmetric states 6 = 0, the asymptotic 
solution to (11. 5p reads 

Q''\u) = g^°)(«)gL°'(-f ) + QL°'(«)Q^'(t) , (4.5) 

in terms of the solution to the two-term recursion relations (14.41) written with the help of the 
roots 6k of the transfer matrix t^°\u) = 2{u — 5i){u — 62) ■ ■ ■ {u — 6l) [37] , 



L 



Now recall that we are interested only in the trajectory with the lowest energy only. The latter 
does not depend on the twist of the operator, i.e., it is L— independent. The reason for this being 
that for the corresponding state only the quadratic Casimir ^2°* is large while all other integrals 
of motion become anomalously small. For the roots of the transfer matrix this is translated into 
the statement that just two roots Si = 6l are much larger than the rest of 5's which are negligible 
[57] . yielding the relation 

Sf^-qf'/2. (4.7) 

In this case the genus— (L — 2) hyperelliptic Riemann surface parametrizing the magnet, with 
its moduli determined by the conserved charges q'j^\ degenerates into a sphere, i.e., the spectral 
curve of twist-two operators [37] . This implies that all zones but one of allowed classical motion 
in separated variables collapse into points. In this limit the transfer matrix reduces to t^'^^{u) ~ 
u^T^'^^{u) = u^{2 — N'^/u'^) and the solutions to the recursion relations (14. 4p becomes symmetric 
under the interchange u —>■ —u and equal, Q^^{u) = Q^^\u). In the infinite-spin limit, we then 
find that the leading behavior of the Baxter function is 

{i In Qf{u))' = ij{-iu + i6i) + tpi-iu - iSi) + ...^21nAr+..., (4.8) 

where in the last step we imposed the condition that the evaluation of the anomalous dimensions 
(11.41) requires u ~ and thus it can be neglected compared to A^. This consideration immedi- 
ately suggests that for the minimal-energy trajectory in the single-logarithmic asymptotics the 
dressing factors in the left-hand side of Eq. (14.41) are irrelevant. Thus they can be reduced 
to — *• and cancelled with the factor extracted from the transfer matrix t^'^\u), making the 
equation L— independent, as expected. The latter is clearly seen in the quasiclassical approach 
when one assumes the spectral parameter to scale with N, i.e., u = Nu and m ~ 1. We will use 
the same argument below to write the all-loop Baxter equation for the lowest trajectory. 



4.2 Beyond one loop 

Let us find the equation for the minimal trajectory starting from the multi-loop Baxter equation 
(12.31) . Again, we have to separate only terms which generate leading behavior in the large- 
spin limit. The transfer matrix degenerates on the minimal trajectory to the one of twist-two 
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operators, i.e., t{u) u^^'^(2u'^ + qiu + ^2)- Notice however that only 0{g^) contributions to the 
charges qi,2{g) can induce the leading effect in the large— limit since the quantum corrections 
grow at most logarithmically with 00. Therefore, we can replace t{u) ^ t^°''{u) in the 

right-hand side of (12. 3p . Hence the reduced Baxter equation admits the form 

e"+(^+)Q(M + e"-(^-)Q(M - i) = t^"\u)Q{u) . (4.9) 

Introducing again the ratio of the Baxter functions Q analogous to Eq. (14. 2p . we can write 
again two asymptotic equations for the two components of Q. However, since we are interested 
solely in the lowest trajectory, both equations generate the same contributions to the anomalous 
dimension. Therefore, we may consider only one of the resulting equations, e.g., 

e"+("+)g(M + i)= T^°\u)Q{u) . (4.10) 

Next, introducing the one- and all-loop Hamilton- Jacobi functions, 

S^''\u) = lng'°)(M) , S{u) = InQiu) , (4.11) 

Eq. (14.101) can be rewritten by virtue of the one-loop degenerate Baxter equation (14. 4p for the 
lowest trajectory as follows 

S{u + S^'^u + S{u) + S^'\u) + (T+(x+) = 27rim . (4.12) 

Here m displays the ambiguity in choosing the branch of the logarithm. Since the anomalous 
dimension (13. 5p is expressed in terms of the derivative of the Hamilton- Jacobi function, it is 
instructive to differentiate both side of Eq. (14.120 with respect to u. Using the perturbative 
decomposition of the all-order Hamilton- Jacobi function 

00 

S{u) = S^'\u)+g'S,iu) , with S,{u) = J2g'^-~''>S^-\u) , (4.13) 

n=l 

and rescaling S*^ by extracting its single logarithmic behavior 

iS',,{u) =J:{u)\nN , (4.14) 
we finally arrive at the equation for the cusp anomaly 

S(u + ^) - S(«) + -r^= f [2 + /S(t - gt)] = . (4.15) 

^yui - g^ J -I IT u+ + gt 

The cusp anomalous dimension is then found in terms of S making use of Eq. (13.50 as 

T,^M=9' + 9'j^-Vr~t^J:{^,-gt) . (4.16) 

As we will demonstrate below, there exists yet another expression for the cusp anomalous dimen- 
sion in terms of the rescaled Hamilton- Jacobi function E which leads to realization of an iterative 
perturbative structure of Fcusp in gauge theory. Eqs. (I4.15P and (I4.16P are the main results of 
this sections. If one shifts the spectral parameter as m — *• m — |, one immediately realizes that the 
first two terms give the imaginary part of E for real u. Then the use of a dispersion relation for 
the rescaled Hamilton- Jacobi function in the last term allows us to bring the equation into the 
form of a Fredholm equation of the second kind. Then the large— x asymptotics of the solution to 
this integral equation yields the cusp anomaly 2x[u] QmSh{u + ^)\xiu]^oo = — [rcusp(5')/5'^] In A^. 
However below we choose a slightly different route to solve Eq. (14.150 at weak coupling. 
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5 Weak-coupling expansion 



We will seek the solution to the cusp equation (14.151) in the form [2], [20 



du w*"-^ cu-™-i S ( In - 



(5.1) 



with u = 1 — u. This integral representation immediately diagonalizes the difference terms. The 
change of variables to p = Inu/u brings Eq. (15.11) into the form of a Fourier transform. However, 
before we proceed with the above transformation of Eq. (I4.15p . we will manipulate it at first. We 
notice that the first term in the infinite series expansion in Chebyshev polynomials in Eq. (13. ip is 
determined by the all-order anomalous dimension. Therefore, we can separate it from the kernel 
and rewrite the equation for the cusp anomaly Fcusp in a form which immediately suggests yet 
another relation of the Hamilton- Jacobi function to the cusp anomalous dimension. Performing 
these steps, we find 



sinh ) E(p) + 



1 — Ji{,9P) + IT 



dp'- -\]{gp,gp')T.{p')=Q 



V/2 



P — p 



(5.2) 



where the kernel U is expressed in terms of the Bessel functions. 



vip,p') = Mp) 



Mp')--Mp') 
p 



Mp) - -Mp) 
p 



(5.3) 



An examination of Eq. (15. 2p immediately suggests that the last term dies out for p — much 
faster than the first two, which scale linearly with p. Therefore, we deduce yet another represen- 
tation for Fcusp in terms of the solution S to the cusp equation (15. 2p . namely. 



E(0) 



cusp 



9' 



(5.4) 



At the same time, we can use Eq. (14.161) for the anomalous dimension in terms of the Hamilton- 
Jacobi function, such that we get 
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-e~^/Vi((7p)S(p). 
P 



(5.5) 



This expression clearly displays the mixing of orders and thus exhibits an iterative structure 
of the perturbative series in coupling constant, i.e., the cusp anomaly at higher orders can be 
determined in terms of S(p) at lower orders. Combining Eqs. (15.21) - (15.51) together we reproduce 
the cusp equation derived in Ref. [T2\ . 

Finally, let us solve the cusp equation perturbatively. Writing the expansion in coupling 
constant as 

sraf:^^"^"(-)' 



(5.6) 



n=0 



where the prefactor is extracted for the latter convenience, and substituting it into the cusp 



8 



equation l\5.2h . we find for the few lowest order functions 



So(p) = -1 , 

MP) = - ^ - f4c(5) + ^co)! p + - Imp' + ^ + 



20160 8 V16 96'^ 7' 480' 96'^'^ 2304 9216 

Tlie p-independent term in tliese expressions determines tlie cusp anomaly according to Eq. (15 ■4p . 
The lowest six orders of Fcusp read 

rcusp(^7) = 9'- ^g' + ^4^9' (5.7) 



12" 720 



TT 



2 



20160 8 r \ 907200 48 



C(3)^-oC(3)C(5)b 
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/1368837ri° ttV^,, Stt^ 51 2 , ^2 



- — C(3)' - — C(3)C(5) - — C(5)' - — C(3)C(7) ^ 

V479001600 240^^ 43^^^^^^ 64^^^ 64^^ ^^^7^ 

7680089." _ ((3)^ _ ll£!^(3)^(5) _ U^^^.y- _ ^c(3)C(7) 

V87178291200 48384^^^ 64 1920^^^^ 128^^ 128^^^^ 

^^^-C(5)C(7)-^C(3)C(9)V'' + ---- 



The two- and three-loop coefficients agree with Feynman diagram calculations of Refs. |22l 
and [251 EH Us] , respectively, and the rest with available predictions of Ref . [12j. The calculation 
can be extended to few dozens of terms in the series fl5.6p . but the results are too cumbersome 
to display here. 



6 Outlook 



In this note we proposed a multi-loop asymptotic Baxter equation for anomalous dimensions of 
arbitrary twist— L, spin— single-trace holomorphic Wilson operators in maximally supersym- 
metric Yang-Mills theory. We developed an approach for the asymptotic solution of the resulting 
equation for large values of spin and derived an all-order equation for the cusp anomaly which 
governs the Sudakov asymptotics of anomalous dimensions. The problem with the asymptotic 
nature of the equation was overcome by studying the lowest-energy trajectory which is insensitive 
to the twist of the operator in the single logarithmic regime Le^ -C iV, L, iV — > oo. 

There are many questions which remain to be addressed. One has to constrain the amount 
of ambiguity left in restoration of higher loop effects from the lowest few terms of perturbative 
series for the dressing factors. The analysis of the strong-coupling expansion of Fcusp is of special 
interest in light of available predictions for it from string theory [3l]. A preliminary analysis 
reveals however that (7 = cxd is an essential singularity of the cusp equation. Next, one has to 
understand how to incorporate wrapping effects to the Baxter equation (12.31) and to identify a 
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putative microscopic spin chain standing behind it. An ultimate goal would be to generalize the 
all-order Baxter equation to all sectors of A/" = 4 super- Yang-Mills theory which is conceivably 
described by a long-range graded magnet. 



Note added 

Recently a new calculation was published of the four-loop cusp anomalous dimension using the 
unitarity technique [12]. Their numerical finding explicitly demonstrates that the prediction (15.71) 
based on the Baxter equation (12. 3p with the dressing factor (13.11) is incorrect starting from four 
loops. In a companion paper [13], a modified form of the cusp equation was proposed which 
takes into account a non-trivial dressing factor in Bethe equations of Ref . [2l] . 

Presently we use the result of Ref. [^ H3] in order to fix the form of the four-loop correction 
to the Baxter equation (12. 3p and find anomalous dimensions of local Wilson operators. It was 
suggested [12], that to reconcile within error bars the result of their numerical calculation with 
the one coming from the cusp equation, the sign of the C^(3) in four-loop contribution of Eq. 
(15. 7p has to be flipped. This requires the following additive modification of the four-loop cusp 
anomaly (15. 7p . 

T,.sA9) = -^9' + Oig''). (6.1) 



In order to generate it from the cusp equation, one has to add the following terno to the left-hand 
side of Eq. f lO) 

■■■ + 2ag'Mgp) {l + g dp' e'^' l't{p')^-^^ + 0{g') , (6.2) 

in agreement with Ref. [21]. Here the favored value of the constant is a = |C(3) 112] • This 
translates into a modification of the integrand in Eq. (I4.15p . 

' ' + (6.3) 



u+ + gt u+ + gt 

A simple analysis allows to unambiguously restore the correction term to the dressing factors 
(13. ip of the Baxter equation (12.30 . Namely, the former get shifted as 

(T± (x) (T± (x) + A± (x) , (6.4) 

with 

A±(x) = ^ VT^ [lnQ(±f - gt)]' + 0{g') . (6.5) 

Taking into account this extra term, the anomalous dimensions of Wilson operators acquire 
additional contributions. For instance, the four-loop term in Eq. (13.60 gets corrected by 

l{9) = ■ ■ ■ - a^-^g' + 0{g'') . (6.6) 



"'in the unnumbered equations above ()5.2p . it yields corrections to the right-hand side the equations, i.e., 
S2(p) = • • • + ^ap, E3(p) = • • • + ^aC(3) + ^ap(7r2 
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This explicitly demonstrates that the attempt to rescue the principle of maximal transcendental- 
ity [26] in the cusp anomalous dimension with a = |C(3) results in breaking of the rational form of 
anomalous dimensions of local Wilson operators, i.e., they acquire transcendental addenda (16.61) 
in addition to rational terms (13.61) . At the current state-of-the-art of higher loop calculations 
such terms are not ruled out yet. Within Mueller's cut vertex technique the main sources 
of transcendental constants in local anomalous dimensions comes from virtual self-energy and 
vertex corrections with rational terms being generated by real cuts. The finiteness of maximally 
supersymmetric Yang-Mills theory, especially transparent in the light-cone gauge where Ward 
identities imply equality of the vanishing beta function with all field renormalization constants, 
seems to suggests the absence of transcendental constants in local anomalous dimensions. This 
question deserves however a thorough study. 



This work was supported by the U.S. National Science Foundation under grant no. PHY-0456520. 
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